Chaotic Spin Correlations in Frustrated Ising Hierarchical Lattices 
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Spin-spin correlations are calculated in frustrated hierarchical Ising models that exhibit chaotic 
renormalization-group behavior. The spin-spin correlations, as a function of distance, behave chaot- 
ically. The far correlations, but not the near correlations, are sensitive to small changes in tem- 
perature or frustration, with temperature changes having a larger effect. On the other hand, the 
calculated free energy, internal energy, and entropy are smooth functions of temperature. The 
recursion-matrix calculation of thermodynamic densities in a chaotic band is demonstrated. The 
leading Lyapunov exponents is calculated as a function of frustration. 
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It was shown some time ago that frustrated Ising spin 
magnetic systems exhibit chaotic behavior of the in- 
teraction constants under renormalization-group trans- 
formations, which readily leads to the description of a 
spin-glass phase. [l| This chaotic rescaling behavior was 
originally demonstrated in frustrated, but non-random 
systems. It was subsequently shown that the same 
chaotic rescaling behavior occurs in quenched random 
spin glasses. @ Chaotic rescaling behavior has now been 
established as the signature of a spin-glass phase. [EEL IE 
IE SB BEE III 111 Although the chaotic behavior of the 
interaction constants was demonstrated in frustrated sys- 
tems, the behavior of spin-spin correlation functions and 
the instabilities to initial conditions had not been calcu- 
lated to-date. This study presents such results, yielding 
both smooth and unsmooth behaviors, as seen below. In 
this process, the recursion-matrix calculation of thermo- 
dynamic densities in a chaotic band is demonstrated. 




FIG. 1: The family of hierarchical models from Ref.[l|. In 
unit (a), there are p cross bonds. In unit (b), two paths, 
consisting of mi and rri2 > mi bonds in series, are in parallel. 
In (c), the final graph of the model is assembled with p a and 
Pb of each unit in parallel. Each wiggly bond, representing an 
infinite antiferromagnetic coupling, has the effect of reversing 
the sign of an adjoining bond. 



Hierarchical models are exactly soluble models that ex- 
hibit non-trivial cooperative and phase transition behav- 
iors [IE EE 13 and have therefore become the testing 
grounds for a large variety of phenomena, as also seen in 
recent works. [H EE EE OS Efl EH, El, HE HE EE P . 
The hierarchical models in which the chaotic rescaling 
behavior of the interaction constants under frustration 
was seen [lj are defined in Fig[T] The two units [FigQJa) 
and (b)] assembled in the construction of these lattices 
a priori represent the generically distinct local effects 
of frustration occurring in spin-glass systems on conven- 
tional lattices: In FigfTJa), correlation at the small scale 
(vertical bonds) inhibit at low temperatures the propaga- 
tion of correlation at the larger scale (horizontally across 
the unit), namely causing a disordering by ordering. In 
Fignjb) , competition between paths of different lengths 
weakens but does not eliminate the propagation of cor- 
relation across the unit. These two generic effects are 
incorporated into the hierarchical lattices of FigfT] No 
other such generic effects occur in spin glasses. 



Renormalization- Group Transformation. Hierarchical 
lattices are constructed by the repeated self-imbedding 
of graphs. fl3l. ITU Their solution, by renormalization- 
group theory, consists of the reverse procedure. The 
number of bonds of the imbedding graph gives the vol- 
ume rescaling factor, b d = (4 + p)p a + (mi + m2)pb in 
the current case, and the shortest path length across the 
imbedding graph gives the length rescaling factor, 6 = 2 
here, leading here to a dimensionality d that is greater 
than 2. Each straight line segment in FigQ] corresponds 
to an interaction —/3Hij — KoiO-j + G with K > be- 
tween Ising spins aj = ±1 at vertices i. Frustration is 
introduced by the wiggly bonds. The additive constant 
G is generated by the renormalization-group transforma- 
tion and enters the calculation of the thermodynamic 
functions and correlations of the original, unrenormalized 
system. [29(. The renormalization-group transformation 
consists in summing, in the partition function, over the 
internal spins of the innermost imbedding graphs, which 
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FIG. 3: (a) Number of visits per interaction interval At — 
0.005, for 5000 chaotic iterations in the trajectory starting at 
t (0) = 0.5, for p = 4, p a = 40, pi, = l,m\ = 4.7, m 2 = m 1 + 1. 
(b) Number of visits per correlation interval A < SiSj > = 
0.005 for the trajectory in (a) 



FIG. 2: Renormalization-group flow topologies, Lyapunov ex- 
ponents, and spin-spin correlations, as Hamiltonian parame- 
ters are scanned. In both examples, p = 5,pb ~ 2, mi = 
mi + 5. In the upper panels, the lower curve, visible on the 
left, is a line of unstable fixed points, giving the second-order 
phase transition between the paramagnetic (below) and or- 
dered (above) phases. In the ordered phase, only fixed points, 
limit cycles, and chaotic bands that are stable (attractive) are 
shown. In the lower panels, the spin-spin correlation function 
< SiSj > for spins separated by a distance 2™ are given for 
consecutive n at each value of the Hamiltonian parameters. 
It is shown in the text that the Lyapunov exponents, middle 
panels, apply to both upper and lower panels. Left panels: 
Scanning p a , which increases the disordering by ordering ef- 
fect, at mi — 5. Right panels: Scanning mi, which increases 
the ground-state entropy of per first renormalized bond, at 
Pa = 50. 



are thereby replaced by a renormalized bond with 

K 1 = Pa tanh" 1 t a + p b (tanh" 1 t mi - tanh" 1 t™ 2 ), 
where t a = 2t 2 {l - f)/(l + t 4 - 2t 2 t), 
G' = b d G + [2 Pa + (mi + m 2 - 2)p b ] In 2 - PaP K 



Pa (1 + t 2 ) 2 - AtH 

2 n (l-ta)2(i_t)2 



p b ln (l-t 2mi )(l-i 2m2 ) 
2 n (1 - t 2 ) mi+m2 



(1) 



with t = tanh if and t = tanh(pif). 

Results: Self-Similar Chaotic Bands. When this fam- 
ily of models is scanned as a function of p or m\, re- 
spectively increasing the disordering by ordering effect or 
the ground-state entropy per first renormalized bond, the 
chaotic behavior of the renormalization-group is entered, 
in the low-temperature phase, via the series of period- 
doubling bifurcations, with the Feigenbaum exponent of 
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tions for the trajectory in Fig. 3(a). The relative difference 

200 " 



5gnSf£i|Ani|/nj, between two consecutive groups, in the 



number of visits 7ij to each of the 200 interaction intervals 
i is shown as a function of group size N. 
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FIG. 5: The spin-spin correlation function < SiSj > for spins 
separated by a distance 2 n , for K = 2.5, p = 4,p a = 40, 
I . mi = 3.7, 4.7, 5.7, mi = mi + 1. 



4.669 as shown in FigJ^l An example of the chaotic 
bands of the interaction constant is shown in FigJ3][a). 
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FIG. 6: Deviations, for small temperature or frustration 
change, in the spin-spin correlation function < SiSj > for 
spins separated by a distance 2 n , for K = 0.8, p — 4,p a = 
40, ph = 1, mi — 8, m.2 = mi + 1. In (a) and (b), between the 
two trajectories, AK = 0.001 and Ap a = 0.001 respectively. 

Discovery of these chaotic bands immediately led to a 
spin-glass interpretation: Under repeated scale changes, 
the entire band is visited by the effective coupling of the 
length scales that are reached after each renormalization- 
group transformation. This chaotic sequence of hopping 
visits stretches from the strong coupling to the weak cou- 
pling edges of the band. This signifies that, as the system 
is viewed at successive length scales, strong and weak 
correlations are encountered in a frozen but chaotic se- 
quence, meaning a spin-glass phase. This interpretation 
had not been followed by an actual calculation of these 
spin-spin correlations, which is done in the current study. 

Also shown in Figj2] are the leading Lyapunov expo- 
nents (A), used to describe the behavior of a dynami- 
cal system that starts at xq and evolves for n iterations, 

Xi+l = f(Xi), 

A = Um IVln|^±i|. (2) 
n -> oo n f-^ dx k 

fc=0 

The iteration function / may depend on different param- 
eters, as our iteration function K'(K) depends on mi and 
p a - Such an iterated map function has a chaotic trajec- 
tory for a particular parameter value if the Lyapunov 
exponent is positive. Conversely, a negative A indicates 
eventual attraction to a fixed point or a limit cycle. A bi- 
furcation point, where a period doubling occurs, is iden- 
tified with A being zero.[28[ Furthermore, 

d < SjSj >k+i 
d < SiSj > fc 

d < SjSj >k+i rftanfi(jffc + i) dtemh(K k ) ^ 

dtaah.(Kk+i) dtanh(Kk) d < SiSj >k ' 
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FIG. 7: For model parameters p — 4,p a = 10, pb = 1, mi = 
7, m,2 = mi + 1, (a) free energy per bond, F — \nZ/Nbond, 
(b) internal energy, U =< SiSj >, (c) entropy per bond, S = 
In Z/Nbond — K < SiSj >, versus temperature K~ . 

so that the first and last factors from Eq. cancel out 
in the successive terms in Eq. Thus, the interaction 
constants and the spin-spin correlations have the same 
Lyapunov exponents. 

An important characteristic of the chaotic bands is 
that they are self-similar: After the transient behavior 
of a number of renormalization-group transformations, 
the profile of the chaotic band formed by each successive 
group of N renormalization-group calculations becomes 
identical in the limit of large N. The overlaps between 
such successively formed bands is shown as a function of 
N in Fig. 4. Physically, this signifies that a geometrically 
coarse-grained spin-glass phase is self-similar. This prop- 
erty of the chaotic bands is important in the calculation 
of the correlation functions. 

Results: Calculation of the Correlation Functions. The 
recursion relations for the densities is [2!| 

/ ud 9CT \ 

[1, < s iSj >} = b- d [l, < S iSj >'] ° OK, . (4 ) 

V U OK ' 

In an ordinary renormalization-group analysis, this den- 
sity recursion relation is iterated, 

[i,<s iSj >] = b- dn [i,< SiSj >(»)]. r* n >-r< n - 1 >-...-r< 1 > ) 

until the (n)th renormalized system is as close as one 
desires to a sink fixed point, and the renormalized den- 
sities [1,< SiSj >( n ^] are inserted as the left eigenvector 
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of TW with eigenvalue 6 d .[29j In the current calculation, 
this cannot be done, since the renormalization-group tra- 
jectory does not approach a sink fixed point, but chaoti- 
cally wanders inside a band. On the other hand, in this 
chaotic-band sink, we obtain the limiting behavior, due 
to the self-similarity property of the chaotic band, 

b -dn T (n) . T (n-1) . _ _ . T (l) / 1 -J ^ ^ (g) 

so that X —< SiSj > and this result is independent of 
the chaotic-band terminus < SjSj >( n >. The disappear- 
ance of the lower diagonal reflects dK'/dK < b d , itself 
due to sequential non-infinite bonds and frustration in 
the chaotic-band sink. Alternately, < s^s,,- > can be cal- 
culated from numerical differentiation of the free energy 
obtained from the renormalization of the additive con- 
stant G. 

The calculated spin-spin correlations as a function of 
spin separation are shown in FigO It is seen that 
the spin-spin correlations behave chaotically, for all dis- 
tances, between strong and weak correlations, numeri- 
cally justifying the spin-glass phase interpretation. Thus, 
spin-spin correlations also span chaotic bands, as illus- 
trated in Fig. 2, lower panels, and Fig. [31(b). 
Results: Unsmooth and smooth behaviors. FigJB] shows 
the behavior, at all distances, of the spin-spin corre- 
lations under small changes in temperature or frustra- 



tion. It is seen that the near correlations are unaffected, 
whereas the far correlations are strongly affected, namely 
randomly changed, with temperature changes having a 
larger such effect. 

Finally, the free energy, calculated from summing 
the additive constants generated by the successive 
renormalization-group transformations, the internal en- 
ergy, calculated from the nearest-neighbor spin-spin cor- 
relation, and the entropy are shown in FigslTJa-c) as a 
function of temperature. They exhibit smooth behav- 
iors. Zero-tcmperature entropy [30[, due to frustration, 
is seen. 

In closing, we note that other forms of chaotic be- 
havior, namely as a function of system size (3lj | or 
as the chaos of near-neighbor correlations in the zero- 
temperature limit for appropriately chosen interactions 
[32| . intriguingly occur in spin-glass systems. In contrast 
to our current results, the renormalization of correlations 
in a strange non-chaotic attractor are given in (33| . 
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